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(Dated: December 13, 2006) 

A measurement of the top quark pair production cross section in proton anti-proton collisions at 
an interaction energy of y/s — 1.96 TeV is presented. This analysis uses 405 pb _1 of data collected 
with the D0 detector at the Fermilab Tevatron Collider. Fully hadronic it decays with final states 
of six or more jets are separated from the multijet background using secondary vertex tagging and 
a neural network. The tt cross section is measured as a t i — 4.5l 2 'g(stat)l}' 4 (syst) ± 0.3(lumi) pb 
for a top quark mass of m t = 175 GeV/c 2 . 



PACS numbers: 13.85.Lg, 13.85.Ni, 14.65. Ha 
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The standard model (SM) predicts that the top quark 
decays primarily into a W boson and a b quark. The mea- 
surement presented here tests the prediction of the SM 
in the dominant decay mode of the tt system: when both 
W bosons decay to quarks, the so-called fully hadronic 
decay channel. This topology occurs in 46% of tt events. 
The theoretical signature for fully hadronic tt events is 
six or more jets originating from the hadronization of the 
six quarks. Of the six jets, two originate from b quark 
decays. Fully hadronic tt events are difficult to identify 
at hadron colliders because the background rate is many 
orders of magnitude larger than that of the tt signal. 

We report a measurement of the production cross- 
section of top quark pairs, <r £ f, using data collected with 
D0 in the fully hadronic channel, that exploits the long 
lifetime of the b hadrons in identifying b jets. To increase 
the sensitivity for tt events, we used a neural network to 
distinguish signal from the overwhelming background of 
multijet production through Quantum Chromodynamic 
processes (QCD). 

The D0 detector [1] has a central tracking system con- 
sisting of a silicon micro strip tracker (SMT) and a central 
fiber tracker (CFT), both located within a 2 T supercon- 
ducting solenoidal magnet, with designs optimized for 
tracking and vertexing at pseudorapidities \rj\ < 3 and 
1 77 1 < 2.5, respectively. Rapidity y and pseudorapidity 
r] are defined as functions of the polar angle and pa- 
rameter p as 2/(0,/?) = |ln[(l + /3cos0)/(l - /3cos0)] 
and 77(0) = ?/(0, 1), where /3 is the ratio of the particle's 
momentum to its energy. The liquid-argon and uranium 
calorimeter has a central section (CC) covering pseudo- 
rapidities 1 77 1 up to ~ 1.1 and two end calorimeters (EC) 
that extend coverage to \rj\ « 4.2, with all three housed 
in separate cryostats. Each calorimeter cryostat con- 
tains a multilayer electromagnetic calorimeter, a finely 
segmented hadronic calorimeter and a third hadronic 
calorimeter that is more coarsely segmented, providing 
both segmentation in depth and in projective towers of 
size 0.1 x 0.1 in rj-c/) space, where (j) is the azimuthal an- 
gle in radians. An outer muon system, covering \rj\ < 2, 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T iron toroids, followed 
by two similar layers after the toroids. The luminosity is 
measured using plastic scintillator arrays placed in front 
of the EC cryostats. 

The data set was collected between 2002 and 2004, 
and corresponds to an integrated luminosity C = 
405 ± 25 pb _1 [2]. To isolate events with six jets, we 
used a dedicated multijet trigger. The requirements on 
the trigger, particularly on jet and trigger tower energy 
thresholds, were tightened during the collection of the 
data set to manage the increasing instantaneous lumi- 
nosities delivered by the Fermilab Tevatron Collider. The 
change in trigger requirements had little effect on the ef- 
ficiency for signal, while removing an increasing number 
of background events [3] . The trigger was tuned for the 
fully hadronic tt channel and was optimized to remain 
as efficient possible while using limited bandwidth. The 



collection rate after all trigger levels was fixed to a few 
Hz, which was completely dominated by QCD multijet 
events as the hadronic tt event production rate is ex- 
pected to be a few events per day. We required three or 
four trigger towers above an energy threshold of 5 GeV 
at the first trigger level, three reconstructed jets with 
transverse energies [Et) above 8 GeV at the second trig- 
ger level, combined with a requirement on the sum of the 
transverse momenta (pr) of the jets, and four or five re- 
constructed jets at transverse energy thresholds between 
10 and 30 GeV at the highest trigger level [l|. 

We simulated tt production using alpgen 1.3 to gen- 
erate the parton-level processes, and pythia 6.2 to 
model hadronization [4, 5]. We used a top quark invari- 
ant mass of m t = 175 GeV/c 2 . The decay of hadrons car- 
rying bottom quarks was modeled using evtgen @ . The 
simulated tt events were processed with the full geant- 
based D0 detector simulation, after which the Monte 
Carlo (MC) events were passed through the same recon- 
struction program as was used for data. The small differ- 
ences between the MC model and the data were corrected 
by matching the properties of the reconstructed objects. 
The residual differences were very small and were cor- 
rected using factors derived from detailed comparisons 
between the MC model and the data for well understood 
SM processes such as the jets in Z boson and QCD dijet 
production. 

In the offline analysis, jets were defined with an iter- 
ative cone algorithm [7]. Before the jet algorithm was 
applied, calorimeter noise was suppressed by removing 
isolated cells whose measured energy was lower than four 
standard deviations above cell pedestal. In the case that 
a cell above this threshold was found to be adjacent to 
one with an energy less than four standard deviations 
above pedestal, the latter was retained if its signal ex- 
ceeded 2.5 standard deviations above pedestal. Cells that 
were reconstructed with negative energies were always re- 
moved. 

The elements for cone jet reconstruction consisted of 
projective towers of calorimeter cells. First, seeds were 
defined using a preclustering algorithm, using calorime- 
ter towers above an energy threshold of 0.5 GeV. The 
cone jet reconstruction, an iterative clustering process 
where the jet axis was required to match the axis of a 
projective cone, was then run using all preclusters above 
1.0 GeV as seeds. As jets from tt production are rel- 
atively narrow due to relatively high jet pr, the jets 
were defi ned using a con e with radius i? C one = 0.5, where 
AR = ^(Ay) 2 + (A(j)) 2 . The resulting jets (proto-jets) 
took into account all energy deposits contained in the jet 
cone. If two proto-jets were within 1 < AR/R cone < 2, 
an additional midpoint clustering was applied, where the 
combination of the two proto-jets was used as a seed 
for a possible additional proto-jet. At this stage, the 
proto-jets that share transverse momentum were exam- 
ined with a splitting and merging algorithm, after which 
each calorimeter tower was assigned to one proto-jet at 
most. The proto-jets were merged if the shared pr ex- 
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ceeded 50% of the pt of the proto-jet with the lowest 
transverse momentum and the towers were added to the 
most energetic proto-jet while the other candidate was re- 
jected. If the proto-jets shared less than half of their pt, 
the shared towers were assigned to the proto-jet which 
was closest in AR space. The collection of stable proto- 
jets remaining was then referred to as the reconstructed 
jets in the event. The minimal pr of a reconstructed jet 
was required to be 8 GeV/c before any energy corrections 
were applied. 

We removed jets caused by electromagnetic particles 
and jets resulting from noise in hadronic sections of the 
calorimeter by requiring that the fraction of the jet en- 
ergy deposited in the calorimeter (EMF) was 0.05 < 
EMF < 0.95 and the fraction of energy in the coarse 
hadronic calorimeter was less than 0.4. Jets formed from 
clusters of calorimeter cells known to be affected by noise 
were also rejected. The remaining noise contribution was 
removed by requiring that the jet also fired the first level 
trigger. 

To correct the calorimeter jet energies back to the level 
of particle jets, a jet energy scale ( JES) correction C JES 
was applied. The same procedure has to be applied to 
Monte Carlo jets to ensure an identical calorimeter re- 
sponse in data and simulation. The particle level or true 
jet energy E true was obtained from the measured jet en- 
ergy E m and the detector pseudorapidity, measured from 
the center of the detector (rjdet), using the relation 

true E m - E (rjdet,C) 



E ] 



7l(ri det ,E™)S(ri deU E™) 



In data and MC the total correction was applied to the 
measured energy E m as a multiplicative factor C JES . 
EoiVdet, £) was the offset energy created by electron- 
ics noise and noise signal caused by the uranium in the 
calorimeter, pile-up energy from previous collisions and 
the additional energy from the underlying physics event. 
The dependence on the luminosity C was caused by the 
fact that the number of additional interactions was de- 
pendent on the instantaneous luminosity, while the de- 
pendence on y was caused by variations in the calorimeter 
occupancy as a function of the jet rapidity. 7Z(r]det, E m ) 
parameterized the energy response of the calorimeter, 
while S(rjdet, E rn ) represents the fraction of the true par- 
tonic jet energy that was deposited inside the jet cone. 
This out-of-cone showering correction depended on the 
energy of the jet and its location in the calorimeter. 

The JES was measured directly using pt conservation 
in photon + jet events. The method was identical for 
data and simulation and used transverse momentum bal- 
ancing between the jet and the photon. As the energy 
scale of the photon was directly and precisely measured 
(the electromagnetic calorimeter response was derived 
from measurements of resonances in the e + e~ spectrum 
like the Z boson), the true jet energy could be derived 
from the difference between the photon and jet energy. 
Eq, 1Z and S were fit as a function of jet rapidity and mea- 
sured energy, which lead to uncertainties coming from the 



C Jhb (E m ,7 ldet ,C)-E m . 

(1) 



fit (statistical) and the method (systematic). The total 
correction C JES was approximately 1.4 for data jets in 
the energy range expected for jets associated with top 
quark events. The uncertainties on C JES ', which were 
dominated by the systematic uncertainty of the out-of- 
cone showering correction Sfjjdet, E m ), were a few per- 
cent and were dependent on the jet energy and rapidity. 

The jet energy resolution was measured in photon + 
jet data for low jet energies and dijet data for higher jet 
energy values. Fits to the transverse energy asymmetry 
[Pt(1) — Pt(2)]/[pt(1) +Pt(2)] between the transverse 
momenta of the back-to-back jets and/or photon (pt(1) 
and Pt(2)) were then used to obtain the jet energy reso- 
lution as a function of jet rapidity and transverse energy. 
The uncertainties on the jet energy resolution were dom- 
inated by limited statistics in the samples used. 

In this analysis, we considered a data set consisting 
of events with four or more reconstructed jets, in which 
the scalar sum of the uncorrected transverse momenta 
fluncorr Q f a jj ^ e jets in the event was greater than 90 
GeV/c. The final analysis sample was a subset of this 
sample, where at least six jets with corrected transverse 
momentum greater than 15 GeV/c and \y\ < 2.5 were 
required. Events with isolated high transverse momen- 
tum electron or muon candidates were vetoed to ensure 
that the all-hadronic and leptonic tt samples were dis- 
joint 11,0]. I n addition, we rejected events where two 
distinct pp interactions with separate primary vertices 
were observed and the jets in the event were not assigned 
to only one of the two primary vertices. The primary 
vertex requirement did not affect minimum bias interac- 
tions or tt events. Table [I] lists the efficiencies after the 
first set of selection cuts, commonly referred to as prese- 
lection, which includes the requirements on the primary 
vertex, the number of reconstructed jets and the presence 
of isolated leptons, and the efficiency after preselection 
and after preselection and the trigger. Besides select- 
ing all hadronic tt events, the analysis was also expected 
to accept a small contribution from the semi-leptonic 
(lepton+jets) tt decay channel. The combined efficiency 
included the fully hadronic and semi-leptonic VF-boson 
branching fractions of 0.4619±0.0048 and 0.4349±0.0027 
respectively [loj . 

We used a secondary vertex tagging algorithm (SVT) 
to identify 6-quark jets. The algorithm was the same as 
used in previously published D0 tt production cross sec- 
tion measurements [1, Q. Secondary vertex candidates 
were reconstructed from two or more tracks in the jet, 
removing vertices consistent with originating from long- 
lived light hadrons as for example Kg and A. Two con- 
figurations of the secondary vertex algorithm were used; 
these were labeled "loose" and "tight" respectively. If a 
reconstructed secondary vertex in the jet had a transverse 
decay length L xy significance (L xy /aL xy ) > 5 (7), the jet 
was tagged as a loose (tight) 6-quark jet. The loose SVT 
was chosen to efficiently identify 6-quark jets, while the 
tight SVT was configured to accept only very few light 
quark jets while sacrificing a small reduction in the em- 
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TABLE I: Efficiency for selection criteria applied before 6-jet 
identification. Efficiencies listed include the efficiency for all 
previous selection criteria. The trigger efficiency is quoted for 
events that have passed the preselection. The uncertainties 
are due to Monte Carlo statistics. Listed are the selection ef- 
ficiencies as determined for tt in the hadronic decay channel, 
the lepton+jets decay channel and the efficiency for all differ- 
ent decay channels corrected for W boson branching fractions. 



cut 



tt — > hadrons tt — > t + jets 



any tt 



preselection 0.2706 ± 0.0016 0.0311 ± 0.0008 0.1385 ± 0.0011 
trigger 0.2527 ± 0.0015 0.0268 ± 0.0007 0.1284 ± 0.0010 



§J1500 




too 1000 % 
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FIG. 1: The Ht distribution for single-tag events (a) and 
double-tag events (b). Shown are the data (points), the back- 
ground (solid line) and the expected tt distribution (filled his- 
togram) multiplied by 140 (60) for the single (double)-tag 
analysis. 



ciency for 6-quark jets. Events with two or more loosely 
tagged jets were called double-tag events. The sample 
that did not contain two loosely tagged jets was inspected 
for events with one tight tag. Events thus isolated were 
labeled single-tag events. The fully exclusive samples 
of single-tag and double-tag events were treated sepa- 
rately because of their different signal-to-background ra- 
tios. The use of the tight SVT selection for single tagged 
events optimized the rejection of mistags, the main back- 
ground in the single-tag analysis. When two tags were 
required, the background sample started to be dominated 
by direct bb production. The choice to use the loose SVT 
optimized the double-tag analysis for signal efficiency in- 
stead of background rejection. 

Compared to light-quark QCD multijet events, tt 
events on average have more jets of higher energy and 
with less boost in the beam direction, resulting in events 
with many central jets that all have similar and rela- 
tively high energies. Moreover, the fully hadronic decay 
makes it possible to reconstruct the W boson and t quark 
four-momenta. To distinguish between signal and back- 
ground, we used the following event characteristics [~Llj : 

(1) Ht'- The scalar sum of the corrected transverse 
momenta of the jets (Fig. 0]). 

(2) E T 6 : The square root of the product of the trans- 
verse momenta of the fifth and sixth leading jet (Fig. [2]). 

(3) A: The aplanarity as calculated from the normal- 
ized momentum tensor (Fig. [3J [1, @, fllj. 
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FIG. 2: The E? 6 distribution for single-tag events (a) and 
double-tag events (b). Shown are the data (points), the back- 
ground (solid line) and the expected tt distribution (filled his- 
togram) multiplied by 140 (60) for the single (double)-tag 
analysis. 
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FIG. 3: The A distribution for single-tag events (a) and 
double-tag events (b). Shown are the data (points), the back- 
ground (solid line) and the expected tt distribution (filled his- 
togram) multiplied by 140 (60) for the single (double)-tag 
analysis. 



(4) [f] 2 ): The ^-weighted mean square of the y of the 
jets in an event (Fig. HJ), see also Ref. [111 ]. 

(5) Ai: The mass-x 2 variable, which was defined as 
M = (M Wl -M w ) 2 /a 2 Mw + (M W2 -M w ) 2 /a 2 Mw + (m tl - 

m t 2 ) 2 / a rn t i wnere the parameters Mw, &m w and a rrif 
were the invariant mass and mass resolution from the jet 





FIG. 4: The (rj 2 ) distribution for single-tag events (a) and 
double-tag events (b). Shown are the data (points), the back- 
ground (solid line) and the expected tt distribution (filled his- 
togram) multiplied by 140 (60) for the single (double)-tag 
analysis. 
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FIG. 5: The M distribution for single-tag events (a) and 
double-tag events (b). Shown are the data (points), the back- 
ground (solid line) and the expected tt distribution (filled his- 
togram) multiplied by 140 (60) for the single (double)-tag 
analysis. 




FIG. 6: The M^ in distribution for single-tag events (a) and 
double-tag events (b). Shown are the data (points), the back- 
ground (solid line) and the expected tt distribution (filled his- 
togram) multiplied by 140 (60) for the single (double)-tag 
analysis. 



four- moment a calculated as observed in all-hadronic tt 
MC, respectively 79, 11 and 21 GeV/c 2 after all correc- 
tions and resolutions were included [12[. Mw { and m ti 
were calculated for every possible permutation of the jets 
in the event. We did not distinguish between tagged and 
untagged jets. The combination of jets that yielded the 
lowest value of Ai is used (Fig. [5|). 

(6) M^ in : The second-smallest dijet mass in the event. 
First, all possible dijet masses were considered and the 
jets that yield the smallest mass were rejected. M^ in 
was the smallest dijet mass as found from the remaining 
jets (Fig. El). 

The top quark production cross section was calcu- 
lated from the output of NN, an artificial neural network 
trained to force its output near 1 for tt events and near — 1 
for QCD multijet events, using the multilayer perceptron 
in the ROOT analysis program [13]. The six parameters 
illustrated in Figs. [TJ6] were used as input for the neu- 
ral net. The very large background-to-signal ratio in the 
untagged data allowed us to use untagged data as back- 
ground input for the training of NN, while tt MC was 
used for the signal. Fig. shows the AW discriminant for 
tt signal and multijet background. Although the distri- 
butions for single- and double-tag events were different 




-1 -0.5 0.5 1 

NN 

FIG. 7: The output discriminant of an artificial neural net- 
work (NN) with six input nodes. All distributions are nor- 
malized to area. NN is optimized to distinguish between fully 
hadronic tt Monte Carlo events (signal) and the background 
from multijet production (background) as predicted by the 
tag rate functions. 



due to increased heavy flavor content in the double-tag 
sample, both samples showed a clear discrimination be- 
tween signal and background. 

The overwhelming background also made it possible 
to use the entire (tagged and untagged) sample to esti- 
mate the background. For the loose and tight SVT, we 
derived a tag rate function (trf — the probability for 
any individual jet to have a secondary vertex tag ) from 
the data with N tag s < 1- The trf was parameterized 
in terms of the pt, <\> and y of the jet and the coordi- 
nate along the beam axis (z) of the primary vertex of 
the event, zpy, in four different Ht bins. To predict the 
number of tagged jets in the event, it was necessary to 
correct for a possible correlation between tagged jets. In 
the single-tag analysis the correlation factor was negligi- 
ble, unlike in the double-tag analysis, where the presence 
of 66+jets events in the sample enhanced the correlation 
correction. We corrected for correlations caused by bb 
background by applying a correlation factor dj, that 
was parameterized as a function of the cone distance be- 
tween the tagged jets, AR. Figure [8] shows the number of 
double-tagged events versus AR as observed in data, and 
the distribution as modeled by the trf with and without 
including dj . We considered significantly different func- 
tional forms for the parameterization of dj an d found 
that the choice of parameterization had little effect on 
the shape of the modeled background distribution. 

The probabilities pi were used to assign a weight, the 
probability that the event could have a given number of 



8 



400 



200 



D0 Run II 405 pb 1 



• double-tagged data 

— TRF prediction, no correlations 

— TRF prediction, incl. correlations 

— TRF prediction, incl. correlations 




2 3 4 5 
AR(tagged jets) 



(0 

4-1 

c 

5 

0) 



400 



200 



• data, N, =1 

' tags 

— background 
■expected tf (6.5 pb) 
■expected s+b 



D0 Run II 405 pb 1 



0.6 



NN 



FIG. 8: The performance of the trf prediction on double- 
tag events (points), without including the correlation factor 
dj (dashed histogram), and including dj for two different 
functional parameterizations (solid histograms). 



FIG. 9: The distribution of the NN output variable for 
single-tag events. Shown are the data (points), back- 
ground (hashed band), signal (filled histogram) and sig- 
nal+background (dashed histogram). The vertical line rep- 
resents the used cut of NN > 0.81. 



tags, to every tagged and untagged event in the sam- 
ple. To ensure the trf prediction was accurate in the 
region of phase space outside the "background" peak of 
the neural network, we used the region — 0.7 < NN < 0.5 
to determine a normalization. In this region of phase 
space, the ti content was negligible. A possible depen- 
dence on ti content was studied by the addition and/or 
subtraction of simulated ti events, as was the variation 
of the interval used for the normalization. Outside the 
background peak, the trf predictions were corrected by: 
SFi = 1.000 ± 0.009 for the single-tag analysis, and 
SF 2 = 0.969 ± 0.014 for the double-tag analysis. The 
errors on the normalization were taken into account as 
a systematic uncertainty on the number of background 
events. 

Both the single-tag and the double-tag analysis were 
expected to be dominated by background, even at large 
values of NN. Figures [9] and [10] show the distribution 
for data (points), the Monte Carlo simulation prediction 
for a t t = 6.5 pb (filled histogram), the background pre- 
diction (line histogram) and the signal+background dis- 
tribution (dashed histogram) 0, Q3] • 

The cross section was calculated from the number of 
ti and background candidates above a cut value of the 
NN discriminant. The cut value was chosen to maxi- 
mize the expected statistical significance s / ^ s + 6, where 
s and b were the number of expected signal and back- 
ground events. The signal and background distributions 
were estimated using the trf prediction and ti Monte 
Carlo events [15]. For both analyses, the expected sta- 
tistical significance was about two standard deviations. 



The optimal cut for the single (double)-tag analysis was 
NN > 0.81 (0.78) shown by a vertical line in Figs. [9] and 
[TUl Table HP gives the observed numbers of events (AT* bs ), 
the background prediction (N£ g ) and the efficiency for 
signal (e t t) that can be used to calculate the ti produc- 
tion cross section via: 



n: 



obs 



■Ni 



bg 



Ml 



(2) 



-TRF) 



where i was "= 1" for the single-tag analysis and "> 2" 
for the double-tag analysis. The number of background 
events is predicted using the trf method. It was likely 
that at values of NN close to unity a certain fraction 
of the sample used to predict the background actually 
consists of tagged or untagged ti events, resulting in an 
increased background prediction. The expected ti con- 
tamination of the background sample was corrected by a 
factor s TRF . In the higher value bins of NN, the contri- 
bution from untagged ti events was significant. s TRF was 
estimated by applying the trf on ti MC, and comparing 
the predicted tagging probability for signal to what was 
expected from background. The size of the Monte Carlo 
sample dominates the uncertainty on s TRF . 

Table HJ lists the systematic uncertainties on the esti- 
mate of the number of background events, the selection 
efficiency and the background contamination. The first 
was uncorrelated between the two analyses, while the lat- 
ter two were correlated as they were derived from the 
same Monte Carlo samples. 

For the single-tag analysis, the systematic uncertainty 
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TABLE II: Overview of observed events, background predic- 
tions and efficiencies. 



symbol 



value 



observed events 
background events 
ti efficiency 
ti contamination 



N; 



T 

obs 

eg 

e TRF 



observed events 
background events 
ti efficiency 
ti contamination 



495 

464.3 ±4.6(syst) 
0.0242t° ; f 8 (syst) 
0.245 ±0.031(syst) 



N- 2 

e^ 2 
°tt 

E^ 2 



439 

400.2t^(syst) 
0.0254l ° ° 6 7 5 o(syst) 
0.194 ± 0.048(syst) 



uncertainty on the normalization and correlation factors 
SF and Cij. For the double-tag analysis, the contribu- 
tion from the uncertainties due to calibration of the b 
quark jet identification efficiency was an additional sys- 
tematic uncertainty on e t f. These uncertainties were de- 
rived by varying the parameterizations used within their 
known uncertainties. 

The single-tag analysis yielded a cross section of 



(7« = 4.1^;g(stat) 



\ + 1.3 



9 (syst) ± 0.3(lumi) pb. 



(3) 



For the double-tag analysis the measured cross section 
was 



&tt 



4.7^2.5 (stat)jlj^(syst) ± 0.3(lumi) pb. 



(4) 
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FIG. 10: The distribution of the NN output variable for 
double-tag events. Shown are the data (points), back- 
ground (hashed band), signal (filled histogram) and sig- 
nal+background (dashed histogram). The vertical line rep- 
resents the used cut of NN > 0.78. 



on the selection efficiency was dominated by the uncer- 
tainty in the jet calibration and identification, which were 
estimated by varying the parameterizations used by one 
standard deviation. The uncertainty on the background 
prediction was dominated by the uncertainty on the trf 
method and the uncertainty on Strf was due to lim- 
ited Monte Carlo statistics. The uncertainty of the trf 
prediction was comprised from the uncertainties coming 
from the fits of the probability density functions at the 
jet level, the statistics of the background sample and the 



As the single-tag and double-tag analysis were mea- 
sured on independent samples, the statistical uncertain- 
ties were uncorrelated. The uncertainties on the selection 
efficiency were completely correlated. Taking all uncer- 
tainties into account, a combined cross section measure- 
ment of 



4.5i?;g(stat)l J^(syst) ± 0.3(lumi) pb (5) 



was obtained, for a top quark mass of m t = 175 GeV/c 2 . 
For a top quark mass of m t = 165 GeV/c 2 , the cross 
section is cr tt ~(165) = 6.2t^(stat)t^(syst) ± 0.4(lumi) 
pb, while for a top quark mass of m t = 185 GeV/c 2 the 
value shifted down to 0^(185) = 4.3^;|(stat)^;o(syst)± 
0.3(lumi) pb. 

In summary, we have measured the ti production cross 
section in pp interactions at = 1.96 TeV in the fully 
hadronic decay channel. We used lifetime 6-tagging and 
an artificial neural network to distinguish ti from back- 
ground. Our measurement yields a value consistent with 
SM predictions and previous measurements. 
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